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Abstract 

The 3+1 formalism of Thorne, Price and Macdonald has been used to derive the linear two-fluid 
equations describing transverse and longitudinal waves propagating in the two-fluid ideal collisionless 
plasmas surrounding a Schwarzschild black hole. The plasma is assumed to be falling in radial direction 
toward the event horizon. The relativistic two-fluid equations have been reformulate, in analogy with 
the special relativistic formulation as explained in an earlier paper, to take account of relativistic 
effects due to the event horizon. Here a WKB approximation is used to derive the local dispersion 
relation for these waves and solved numerically for the wave number k. 
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1 Introduction 

Black holes belong to the most fascinating objects pre- 
dicted by Einstein's theory of gravitation and are still 
mysterious [JJ. Physicists are grappling the theory of 
black holes, while astronomers are searching for real-life 
examples of black holes in the universe [2J. Recent ob- 
servations of nuclei of galaxies have revealed large mass 
concentrations over relatively small volumes [3] |4j (5] H3 
El El [10] . The only plausible explanation for high com- 
pactness of the mass distribution is that galactic nuclei, 
including the Milky Way, contain a massive or a sup- 
permassive (10 6 — 1O 9 M ) black hole. Classical black 
holes do not emit electromagnetic waves and so it is not 
possible to observe them directly. Their presence is in- 
ferred by observing physical processes occurring in the 
plasma situated around their event horizon. Perturbed 
black holes do emit gravitational waves thus making it 
possible to observe them directly. 

Within 3R S (3 Schwarzschild radii) it is possible to 
have plasma [TTJ [T2j [13j [14] . The theory of general rela- 
tivity and its application to the plasma close to the black 
hole horizon have remained esoteric, and little concrete 
astrophysical impact has been felt. In recent year plasma 
equations and general relativity are usually considered 
together. The Coulomb potential of charge particles due 
to coupling is much stronger than the gravitational po- 
tential and are often neglected in the Newtonian approx- 
imation. But the mean gravitational field for certain as- 
tronomical objects like galactic nuclei or black holes may 
be strong and the observation of magnetic fields indi- 
cates that a combination of general relativity and plasma 
physics at least on the level of a fluid description is appro- 
priate. The plasma in the black hole environment may 



act as a fluid and black holes greatly affect the surround- 
ing plasma medium (which is highly magnetized) with 
their enormous gravitational fields. Hence, it is therefore 
of interest to formulate the plasma physics problem in 
context of general relativity. 

A covariant formulation of the theory based on the 
fluid equations of general relativity and Maxwell's equa- 
tions in curved spacetime has so far proved unproductive 
because of the curvature of four-dimensional spacetime 
in the region surrounding a black hole. MacDonald and 
Thorne [TH [TB] have introduced Maxwell's equations in 
3+1 coordinates, which provides a foundation for formu- 
lation of a general relativistic (GR) set of plasma physics 
equations in the strong gravitational field of black hole 
and the "membrane paradigm" [17] is a good example of 
a formalism which has been developed for this purpose. 

The membrane paradigm is mathematically equiva- 
lent to the standard, full general relativistic theory of 
black holes, so far as all physics outside the horizon is 
taken into account in which the black holes event horizon 
are replaced with a surrogate "stretched horizon" viewed 
as a two-dimensional membrane that resides in three- 
dimensional space and evolves in response to driving forces 
from the external universe. But the formulation of all 
physics in this region turns out to be very much simpler 
than it would be using the standard covariant approach 
of general relativity. In the 3+1 formulation, work con- 
nected with black holes has been facilitated by replacing 
the hole's event horizon with a membrane endowed with 
electric charge, electrical conductivity, and finite temper- 
ature and entropy so that the physics outside the event 
horizon turns out to be very much simpler than it would 
be using the standard covariant approach of general rel- 
ativity. 



Exploiting the 3+1 formalism, a lot of works have 
been carried out. The study of plasma wave in the pres- 
ence of strong gravitational fields using the 3+1 approach 
is still in its early stages. Zhang [T5J [TH] has considered 
the care of ideal magneto hydrodynamics waves near a 
Kerr black hole, accreting for the effects of the holes an- 
gular momentum but ignoring the effects due to the black 
hole horizon. Holcomb and Tajima [20], Holcomb [21] . 
and Dettmann et. al. [22) have considered some proper- 
ties of wave propagation in a Friedmann universe. Daniel 
and Tajima [23] studied the physics of high frequency 
electromagnetic waves in a strong Schwarzschild plasma. 

The most interesting plasma phenomena due to in- 
volvement of collective effects arc strongly nonlinear even 
in the Newtonian or special relativistic case. The addi- 
tion of gravity can be expected to cause even more highly 
nonlinear and violent phenomena. Existing plasma the- 
ories and computational codes are based on Newtonian 
physics, perhaps accounting for special relativity. Sakai 
and Kawata (SK) [23] have developed the linearized treat- 
ment of plasma waves using special relativistic formula- 
tion. Connecting this work Buzzi, Hines, and Treumann 
(BHT) [25] [26] have developed general relativistic two- 
fluid plasma waves around Schwarzschild black hole using 
local approximation. Using local approximation near the 
event horizon of different nonrotating black holes some 
research have been made [27] [28] [29] [30] . Recently, in 
Ref. [31] the two-fluid equations for transverse and longi- 
tudinal waves are simplified using action principle devel- 
oped by Heintzmann and Novello [32] and solved using 
the analytical method developed by Mikhailovskii [3"3"] . 
In this paper, the dispersion relations for transverse and 
longitudinal waves are solved numerically using WKB 
approximation. 

The principal objective of the work presented in this 
paper is to make an initial attempt to be the discovery of 
an instability caused by the general relativistic term in 
the local dispersion relations for transverse (electromag- 
netic) and longitudinal (electrostatic) waves using WKB 
approximation. Such an investigation of wave propa- 
gation in a general relativistic two-fluid plasmas near a 
black hole is important for an understanding of plasma 
processes. That is, what happens when the plasma are 
assumed to be infalling onto the black hole. 

In the present paper Sec. [2] summarizes the 3 + 1 
formulation of Schwarzschild black hole Spacetime. The 
3+1 formulism of the two-fluid plasmas expressing con- 
tinuity, the conservation of energy and momentum, and 
Maxwell's equations in Schwarzschild coordinates are pre- 
sented in Sec. [3] The transformation of Schwarzschild 
metric in Rindler coordinates is also described here. By 
assuming the two-fluid plasma falling toward the radial 
direction of the event horizon, the transverse and longi- 
tudinal components of the two-fluid equations are sepa- 
rated by introducing a new complex transverse fields and 
velocities using Cartesian coordinates in Sec. [4] The lin- 
earized two-fluid equations for the transverse and longi- 
tudinal waves are derived by considering a small pertur- 
bation to the fields and fluid parameters in Sec. [5] In 



Sec. [6]we discuss the way in which the unperturbed fields 
and fluid parameters and their derivatives with respect 
to z depend on the surface gravity and the infall velocity 
from the black hole horizon. Thin layers are considered 
to discover the insstability of two fluid plasma wave near 
the horizon, each with its own appropriate mean value of 
the lapse function and the two-fluid equations are sim- 
plifying using WKB approximation in Sec. [7] In Sec. 
[8] we explain the dispersion relations for the transverse 
and longitudinal waves. The numerical method for solv- 
ing the two-fluid equations is discussed in Sec. [9] The 
numerical results for the Alfven, high frequency electro- 
magnetic, and longitudinal waves are also discussed in 
this section. Finally, in Sec. [10] we present our remarks. 
Here, we use units in which G = c = ks = 1. 

2 3+1 Formalism of 
Schwarzschild Spacetime 

As mentioned in the previous section, our work presented 
in this paper is based on the 3+1 formulation of general 
relativity developed by Thorne, Price, and Macdonald 
(TPM) [Tia[T3[n]. Actually the 3 + 1 approach was orig- 
inally developed in 1962 by Arnowitt, Deser, and Misner 
[34] to study the quantization of the gravitational field. 
Since then, their formulation has most been applied in 
studying numerical relativity [35] . TPM extended the 

3 + 1 formalism to include electromagnetism and applied 
it to study electromagnetic effects near the Kerr black 
hole. So their work has opened up many possibilities for 
studying relativistic effects on plasma around black hole 
horizon in the electromagnetic window. The basic con- 
cept behind the 3 + 1 formulation of general relativity is 
to select a preferred set of spacelike hypersurfaces which 
form the level surfaces of a congruence of timelike curves. 
A particular set of these hypersurfaces constitutes a time 
slicing of spacetime. The hypersurfaces considered here 
are of constant universal time t. In the 3 + 1 formulation, 
the Schwarzschild metric is given by 



ds 2 = g^dx^dx" 



, 2M \ , 2 1 

1 \dt 2 + — 

r J (1-2M/V) 



dr 2 +r 2 dn 2 , (1) 



where dtt 2 = d8 2 + sin 2 



The components x^ denote 



spacetime coordinates and fi, v = 0,1,2,3. The hyper- 
surfaces of constant universal time t define an absolute 
three-dimensional space described by the metric 



ds 2 



idx l dx :> 



1 



(1 - 2M/r) 



dr 2 +T 2 dn 2 . (2) 



The indices i, j range over 1, 2, 3 and refer to coordinates 
in absolute space. The curl and divergence operators in 
the 3+1 set of equations are covariant and can be derived 
from the above metric. The purpose of choosing such an 
absolute three-dimensional space is to isolate the effect of 
the gravitational field from the choice of local Cartesian 
coordinates. 



We consider a set of fiducial observers (FIDOs), i.e. 
the observers remaining at rest with respect to this abso- 
lute space, measure their proper time r using clocks that 
they carry with them and make local measurements of 
all physical quantities. Then their all measured quanti- 
ties are defined as FIDO locally measured quantities and 
all rates measured by them are measured using FIDO 
proper time. For the absolute matric space given in Eq. 
((2|), FIDOs use a local Cartesian coordinate system with 
unit basis vectors tangent to the coordinate lines as 



Of = (I — 2M/r) 1 / 2 ^-,e a = ~,e 



dr' 



rdO' 1 



_1 d_ 

rsmO dip 



(3) 



For a spacetime viewpoint rather than a 3 + I split of 
spacetime, the set of orthonormal vectors also includes 
the basis vector for the time coordinate given by 



d 

rf7 



ld_ 

a dt ' 



(4) 



where a is the lapse function (or redshift factor) defined 
by 



a(r) 



dr 
~dt 



1 



1/2 



(5) 



The lapse function a has an important plays the role of 
a gravitational potential and thereby governs the tick- 
ing rates of clocks and redshifts. We can calculate the 
gravitational acceleration felt by a FIDO from the lapse 
function as follows [251 I2S : 



1 M 

a = — Vina = K&fi 



(6) 



of energy and momentum, and Maxwell's equations in 
3 + 1 formalism. These type of fluid components are 
chosen because the dispersion relations that result from 
the following investigation are valid for either one of the 
two fluids since no assumption can be made upon the 
mass, number density, pressure or temperature of the 
fluids. BHT [25J US] have derived the equation of con- 
tinuity and the conservation of energy and momentum 
and Maxwell's equations in 3+1 split of spacetime. The 
derivation of the equations required for all wave distur- 
bances in the plasma including the general relativistic 
effects are detailed here. The reader is referred to TPM 
H3JHSI7] for this material. In the BHT works [2SJH)] , 
the equation of continuity and Maxwell's equations cou- 
pling with the two-fluid having velocity v s , mass m s , 
number density n s , and charge q s for electromagnetic 
fields (E and B) are given by 



d 

— (7 s n s ) + V • (a7 s n s v s ) 
at 



0. 



and 



V B 

V -E 
dB 

~bt 

<9E 
~dt 



0. 

47TCT, 

-V x (aE), 

V x (aB) - 47raJ, 



with the charge and current densities defined as 



E 



J s q s n s , J = y^^ s q s n s v s , 



(8) 

(9) 
(10) 

(11) 

(12) 

(13) 



Equation © shows that far from the black hole event 
horizon i.e., when a — > 1, the gravitational acceleration 
becomes weak and approaches the Newtonian value for 
flat spacetime. However, near the horizon, the gravita- 
tional acceleration approaches infinity as a — > 0. The 
rate of change of any scalar physical quantity or any 
three-dimensional vector or tensor, as measured by a 
FIDO, is defined by the derivative 



D 
Dt 



ld_ 

a dt 



(7) 



v being the velocity of a fluid as measured locally by a 
FIDO. Since all the quantities are measured locally by 
the FIDO, all the vector quantities are neither covariant 
nor contravariant. The derivation of the gravitational 
accelerations for different spacetime matrices are detailed 
in TPM [H US E] and BHT et al [251156], 

3 Two-fluid Equations 
in 3+1 Formalism 

The two-fluid plasmas situated around the black hole 
horizon are considered as ideal fluids having two compo- 
nents either electron-positron or electron-ion to describe 
the two-fluid equations for the continuity, conservation 



where 7 S is the relativistic Lorentz factor and s is 1 for 
electrons and 2 for positrons (or ions). The lapse func- 
tion a signifies the general relativistic effect around a 
Schwarzschild black hole. The equations for the con- 
servation of energy and momentum, as derived by BHT 
[251 , using above Maxwell's equations for electromag- 
netic fields arc given, respectively, for each species s by 



ld_ 

a dt 



= -V • S s + 2a • S s 



and 
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— ?rr ; - , s — £ s a 
a at 
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(aW,). 



(14) 



(15) 



Where the energy density e s , momentum density S s , and 
stress-energy tensor W% for the electromagnetic field 
have been taken as 



— (E 2 

8tt V 



B 



Wi k = — (E 2 
8tt 



S s = — E x B, 

4tt 



B )g J {E J E + B J B 



-in 



(16) 



For a ideal relativistic fluid of species s in three-dimensions, 
the energy density e s , momentum density S s , and stress- 
energy tensor W 3 s k corresponding to the electromagnetic 
field quantities given in Eq. (|16p are 



^■v 2 ), S s 
+ PsHv k - 



(17) 



where P s is the pressure, and e s is the total energy den- 
sity defined by 



m s n s + P s /(~fg - !)• 



(18) 



The gas constant 7 g take the value 4/3 for T — > oo and 
5/3 for T — > 0. The ion temperature profile is closely adi- 
abatic and it approaches 10 12 K near the horizon [36] . 
Far from the (event) horizon electron (positron) tem- 
peratures are essentially equal to the ion temperatures, 
but closer to the horizon the electrons are progressively 
cooled to about 10 8 — 10 9 K by mechanisms like multiple 
Compton scattering and synchrotron radiation. Using 
the conservation of entropy the equation of state can be 
expressed by 

my- 

The full equation of state for a rclativistic fluid, as mea- 
sured in the fluid's rest frame, is as follows [37J ISH] : 



e = m s n s + m a n a 



P s iH^ iy (im s n 3 /P s ) 
l s n s iH^\im s n s /P s ) 



(20) 



where the H% (x) are Hankel functions. Using Eq. (fl~7j) . 
the energy and momentum conservation equations, Eqs. 
(fl4|) and (fT"5|) . can be rewritten for each species s in the 
following form 

a at act 

+ 7 s g s n s E • v s + 2j 2 (e s + P s )a ■ v s = 0, (21) 



and 



1 d 



jt(E s + P s ) ( ~ + v V V) v, tVP» 



j s q s n s (E + v s x B) + v s j s q s n s ~E ■ v s 



ld_ 

a dt 



P s +7 6 2 (e s + P s )[v s (v s -a)-a] 



(22) 



The fluid velocities and fields in the above equations are 
FIDO measured quantities whereas the fluid densities 
and pressures are measured in the fluid rest frame. Al- 
though these equations are valid in a FIDO frame, for 
a = 1 they reduce to the corresponding special relativis- 
tic equations as given by SK [24] and are valid in a frame 
in which both the fluids are at rest. Since the plasma is 
assumed to be falling in radial direction toward the event 
horizon, the transformation from the FIDO frame to the 
commoving (fluid) frame involves a boost velocity, which 
is simply the infall velocity onto the black hole, given by 



vs = (1 



(23) 



so that, the relativistic Lorentz factor takes the form 
7boost = (1 - v 2 ) 1 ' 2 = l/a. 



The linearized two- fluid equations for transverse and lon- 
gitudinal waves derived in the previous section have the 
terms Since a is a function of Schwarzschild ra- 

dial coordinate r, it is impossible to solve the two-fluid 
equations required for these waves analytically, although 
they form the basis of the numerical procedure. Since 
the plasma is assumed to be falling in the horizon, the 
plasma waves varies with the distance from the horizon. 
The Rindler coordinate system, in which space is lo- 
cally Cartesian, provides a good approximation to the 
Schwarzschild metric near the event horizon. In this co- 
ordinates the two-fluid equations can be solved in a re- 
gion very closed to the horizon and cannot be used to 
the case for extremal black holes because there are no 
Rindler coordinates locally near horizons. Near the point 
(r = 2M, 9 = it/2, <fi — 0) the transformation 



x = 2M{6 - n/2),y = 2M0, z = AMa, 



(24) 



converts the Schwarzschild line element given in Eq. ([T]) 
to 

ds 2 = -a 2 dt 2 + [l - a 2 } " 4 dz 2 + [l - a 2 } ~ 2 (dx 2 

+ dy 2 cos 2 — — ), 
y 2M h 

= ~a 2 dt 2 + (dz 2 + dx 2 +dy 2 ){l + o? A^f }> 

(25) 

which is closely approximated by the Rindler geometry 
with a — z/AM so that 



da 

a(z) = zk, — = k, 
oz 



(26) 



where k = 1/ 4M is the surface gravity of Schwarzschild 
black hole. 



4 One-dimensional 
Wave Propagation 

Since the two-fluid plasma are assumed to be falling to- 
ward the horizon due to the strong gravitational field of 
black hole, the waves can be treated as one-dimensional 
and propagating in the radial direction. The transverse 
and longitudinal parts of the two- fluid equations in Schwarzschild 
coordinates can be separated in a form analogous to that 
used by SK, i.e., to begin with, the velocities and fields 
are split into longitudinal and transverse components. 
Let v SXl v sy , and v sy be the velocity components along 
x, y, and z direction, we define a new complex transverse 
fields and velocities by introducing the complex variables 

v sz (z,t) = u 8 (z,t),v s (z,t) = v sx (z,t) +iv sy (z,t), 
B(z,t) = B x (z,t)+iB y (z,t), 

E{z,t)=E x {z,t)+iE y { Zl t), (27) 
and setting 



'''NT 



■Jsxuy — <JsyJ->x — ~^{ v sB* — V*B), 



v sx E y - v sy E x = -{v s E* ~ v*E), (28) 



where the * denotes the complex conjugate, we could 
write the continuity equation, Eq. (|8]l. in the form 



° , s , 9 . . 
— (7 s n s ) + — (a-y s n s u s ) = 0, 



and Poisson's equation, Eq. (1101) . as 



(29) 



(30) 



Adding the component multiplied by i to the com- 
ponent of the Maxwell's equation, Eq. (fTl"]) . the trans- 
verse equation for the newly defined transverse fields and 
velocities takes the form 

1 dB 

a dt 



, d k . 
-i\— + -\E. 

dz a , 



(31) 



In similar fashion, we obtain from the Maxwell's equa- 
tion, Eq. CEU), 

IdE (d k\ , . . . 

— — = i — + - LB - 47re(7 2 n 2 W2 - Ti^i^i). (32) 
a at \az a J 



Differentiating Eq. (|32|) with respect to t and using Eq. 
dSJ) yields 



d 2 d d 2 



47rea— (71272^2 - «i7i«i) 
dt 



(33) 



The longitudinal component of the momentum conser- 
vation equation, Eq. (|2"2"|) . can be separated out by split 
up into its three vector components and the transverse 
component may then be obtained from the and 
components. The resultant longitudinal and transverse 
components of the momentum conservation equation are 
then respectively given by 



Ps-TTT- = ^nsls [E z + ^ (v s B* - v*B) 



Dt 



(1 - u s )p s a - u s q s n s j s ~E ■ v 



1 0P S 
a dt 



dP s 

dz 



(34) 



and 



Dvs 
' Dt 



q s n s j s (E - \v s B z + iu s B) - u s v s p s a 



1 dP s 



v s q s n s j s E ■ v s 



dt 



(35) 



where E • v s = \ {Ev* s + E*v s ) + E z u s and p s is the total 
energy density defined by 



p s = "f 2 (s s + P s ) = 7 2 s {m s n s + T g P s ) 



(36) 



with T g — jg/("f g — 1). In order to investigate the trans- 
verse electromagnetic waves it is more convenient to work 
from a combination of the transverse components of the 
Maxwell's and momentum conservation equations, Eqs. 
(j3"Tj) - ([33"|) and ([35]) . The longitudinal waves can be inves- 
tigated by combining the longitudinal components of the 
equation of continuity, Eq. (|29[) . Poisson equation, Eq. 
(|30[) . and the conservation of momentum equation, Eq. 



5 Linearized Equations for Trans- 
verse and Longitudinal Waves 

We linearize the sets of the two-fluid equations for the 
transverse and longitudinal waves by considering a small 
perturbation. We introduce the quantities as given in 
BHT of the form 

u s (z, t) = u 0s (z) + 5u s (z, t), 

n s (z, t) = n 0s (z) + Sn s (z, t), 

v s (z,t) = Sv s (z,t), P s (z,t) = P 0s (z) + 5P s (z,t), 

Ps(z,t) = p 0s (z) + Sp s (z,t), E(z,i) = 6E(z,t), 

B z (z,i) = B (z) + 8B z (z,t), B{z,t) = 5B{z,t). (37) 

Here, magnetic field has been chosen to lie along the 
radial e z direction. The relativistic Lorentz factor is also 
linearized such that 7 S = 70s + 6j s , where 



70s = (1 - u^) 2 , 5~f s = 7o s u 0s • Su s 



(38) 



Neglecting the product of perturbation terms the conser- 
vation of entropy, Eq. (|1~9")) . is linearized to 



Sn s 



(39) 



Also the total energy density given in Eq. (|3"6"|) is lin- 
earized to 



8p s 



n 0s 



POs 



Sn s +2u 0s ^ s p 0s Su Sl (40) 



where p 0s = ^ s (m s n Qs + T g P Qs ). 

The transverse part of the Maxwell's equations, Eqs. 
(|3"Tj) and (1521 . and their resultant equation, Eq (l3U|) . are 
linearized to give 

IdSB .(d k\ 

ET = - 1 [ez + ~ ) SE > 41 

a dt \dz a J 

IdSE .(d k\. _ , , 

1 1 — H — )oB — 47re(72n 2 U2 - 7i"i^i), 



dt 
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(42) 



d d 2 



ol—^ + 3an- — + k 2 ) SE 



dz 



dz dt 2 



iirea [ n 02 702 — no 1701-57- ■ ( 43 ) 



dt 



dt 



Linearizing the transverse part of momentum conserva- 
tion equation, Eq. (|35|) . differentiating it with respect to 
t and than substituting for the magnetic field using Eq. 
|), we obtain 



d d 
auos^r- + 77 - u 0s K 
dz dt 



iaq s ^osn sB \ d8v 



POs 



dt 



au 0s — + — + u 0s k ) 5E = 0. (44) 

Pos \ dz dt 



The continuity Eq. (|29|). and Poisson's Eq. ((30]) arc 
linearized to obtain 



d 



8 



wq 4 



d_ 

di 



d ( 1 drios 
+ a- — h k + a ; — 

Oz \no s dz 



du 0i 
dz 



Su s = 0, 
(45) 



and 

dSE z 
dz 



47re(n 02 702 - n i7oi) + 47re(7o 2 <5n2 - 7oi#"i) 



47re(n 02 u 27o2' 5 '"2 - n iMoi7oVwi)- 



(46) 



In similar fashion, the longitudinal part of the momen- 
tum conservation equation, Eq. (|34j) . is linearized to give 




0. (47) 



6 Dependence of the Unperturbed 
Values on z 

Since the plasma is falling toward the event horizon, the 
unperturbed radial velocity for each fluid species as mea- 
sured by a FIDO along is assumed to be the infall 
velocity given by 

uos(z) = v s (z) = [l-a 2 (z)]i. (48) 

From the continuity equation, Eq. (|29|) . it follows that 



r 2 aj 0s n Qs u 0s = const. = r 2 H a H lHn H u H , 

where the values with a subscript H are the limiting 
values at the event horizon. The infall velocity at the 
horizon becomes unity so that uh = 1. Since u 0s = vg, 
70s = 1/a; and hence a~fos — &hJh — 1- Also, because 
Vfi = (rjj/r)?, the number density for each species can 
be written as follows: 



n 0s (z) = n Hs Vg(z). 



(49) 



The equation of state, Eq. (|13p. leads to write the un- 
perturbed pressure, 



pt>.(z) = p H .i—y>, 



(50) 



which in terms of the infall velocity can be written as 

Pos(z) =P Hs vp g (z). (51) 

Since Pq s = kgn^gT^g, then with kg = 1, the tempera- 
ture profile is 



T 0s = T Hs vl^- 1] {z). 



(52) 



The unperturbed magnetic field is purely in the radial z 
direction. It does not experience effects of spatial curva- 
ture. From the flux conservation equation V • Bo = 0, it 
follows that 

r 2 Bo(r) = const. 

from which one obtains the unperturbed magnetic field 
in terms of the infall velocity to the form 



B (z) = B H v^{z). 



(53) 



The acceleration of the two-fluid toward the horizon can 
be determined directly from Eq. (|48l) which approxi- 
mated by the Rindler geometry as 



dvfi 
dz 



1 



(54) 



Similarly, using Rindler geometry the derivatives of the 
unperturbed fields and fluid quantities with respect to z 
in terms of infall velocity become 
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dz 
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The equation of continuity, Eq. (|45p . and momentum 
conservation equation, Eq. (|47|) , with the help of Eq. 
(l55l) in the Rindler coordinates become 
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and 
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where v^ a = 2jgj^ s Po s / pg s is the thermal velocity of 
the fluids of species s. Since the 7q s factor involved in 
the energy density pos in the denominator, the 7q s factor 
in the numerator cancels out and therefore, the thermal 
velocity of the two- fluid plasmas are frame independent. 



7 The WKB Approximation 



and 



The work presented in this paper is devoted to the deriva- 
tion of the local dispersion relations for the transverse 
and longitudinal waves propagating in two-fluid plasmas 
falling in radial direction toward the event horizon of a 
Schwarzschild black hole. The main result of the present 
work seems to be the discovery of an instability caused 
by the general relativistic terms in the dispersion rela- 
tion as shown by Sakai and Kawata [24] for the electron- 
positron plasma in the special relativistic case. Due to 
the Einstein's fundamental principle of general relativity, 
all physical process confined to a small laboratory (i.e., 
localized processes) , we here consider the local values of 
lapse function a ~ ag in the range < a Q < 1 so that 
the general relativistic effects can be shown by build up 
a more complete picture by considering a large number 
of layers within this range of ao values. In this approx- 
imation, the lapse function takes ao — at the event 
horizon and ao = 1 for special relativity. 

Free falling in any gravitational field occur exactly 
in the same way as they do in a laboratory outside of 
any gravitational field. Therefore, if no local instability 
is present in special relativistic plasmas, then no local 
instability should occur in a frame of reference comov- 
ing with a plasma free falling into a black hole. Since 
the transformation from this comoving frame to a coor- 
dinate frame of the metric given in Eq. ([T]) is a simple 
Lorenz boost with the velocity uq s — Vff in the radial 
direction, no local instability should also appear without 
considering non-local effect. Perturbations can grow or 
decay due to non-local effects related to the general rel- 
ativity as well as to the fact that the initial equilibrium 
state of the plasma is non-uniform. We assume that the 
wavelength is small compared with the range over which 
the equilibrium quantities change significantly. Then the 
wavelength must be smaller in magnitude in comparison 
with the scale of the gradient of the lapse function a. 
The mathematically correct way to include non-local ef- 
fects is to consider WKB approximation, that is to write 
all variables as V = V a (z)exp(i J k(z)dz — itot), where 
V a {z) and k(z) are functions, which vary slowly on the 
scale of the wavelength A = 2ir/k(z), i.e. Xdk/dz <C 1 
and XdVa/dz < 1. 

Since the only scale in this work is the radius of the 
black hole horizon, m- Taking the Fourier transforma- 
tion, the set of two-fluid equations for transverse electro- 
magnetic waves, i.e., Maxwell's Eqs. (|41[) -(j43 |) and the 
transverse part of the momentum conservation Eq. (|44l) 
reduced to give 

(a k(z) + in) SE + iuiSB = 0, (58) 

iujSE = (a k(z) + in) SB + 4naoe (702 no2#«2 — 701^01^1), 

(59) 

(a^k{z) 2 — i3aonk(z) — uj 2 — k 2 ) SE 

= i47rea w(n 2702<5w2 - rioi7oi^«i), (60) 
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(aok(z)uo s — lu — iuQ S K) SE = 0. (61) 



Similarly, the set of two-fluid equations for longitudi- 
nal waves described by Eqs. (|46|) . (f56|) . and (|57|) when 
Fourier transformed, become 
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and 



ik5E z = 47re(n 2702 - noiToi) + 47re(7 2<5ri2 - 7oi<5m) 

+47re(n 2Wo27oV u 2 - n iuoi7oi^ u i)- ( 64 ) 

Here, the terms arising from the derivatives dV a /dz and 
dk/dz in the above two- fluid equations for transverse and 
longitudinal waves have been neglected. As in any stan- 
dard WKB approach, the ratio of these neglected terms 
to the main term is ~ A/rjf ~ l/(krn) <C 1. Since 
we have already neglected this terms, it is not mathe- 
matically correct to keep any other terms of the order 
of - X/r H in Eqs. ([58 > (P ]) . Since k - l/{2r H ) < 1, 
we need to neglect the terms containing k compared to 
the main term a$k. Thus, the set of transverse two-fluid 
equations, Eqs. (|58 l> - ([6T1) become 
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iloSE = a kSB + 47ra e (702^02 Sv 2 - 701^01^1), (66) 
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Similarly, the longitudinal two-fluid Eqs. (R)2"j) and (|63p 
follows, as above, to become 

n s7os ( a ok - u 0s uj) 5u s + (a u 0s k - uj) Sn s = 0, (69) 

and 
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8 Dispersion Relations for Trans- 
verse & Longitudinal Waves 

The dispersion relation describing Alfven and high fre- 
quency electromagnetic waves propagating parallel to the 
unperturbed magnetic field Bq may be obtained by elimi- 
nating 5v\, 5i)2 , and SE from the transverse parts of Eqs. 
(1671) and JBSJ as 
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(71) 



for either the electron-positron or electron-ion plasma. 
Where the local plasma frequency for transverse wave 

is Wp 8 = ^/ 47re J° s ~ njL i which depends on the local num- 
ber density of electrons and also the local value of the 
lapse function ao, and the local cyclotron frequency is 

2 d 

oj cs = e i°i™°s a o ^ w \yi C \y depends upon both the local 
magnetic field and the lapse function. Note that, the 
plasma frequency and cyclotron frequency are frame in- 
dependent (independent of 7 0s ) as the thermal veloc- 
ity of the fluids. This is because a boost to the fluid 
frame involves the transformation Bq — > 70s -Bo and po s — 
7g s (m s n 0s + TgP 0s ). We see that general relativistic ef- 
fects enter this dispersion relation only in the ratio cj/ao- 
If one uses local time dr = agdt of the falling observer 
instead of a global coordinate time t, then the dispersion 
relation is reduced to the special relativistic version as it 
should be according to the Einstein relativity principle. 

Similarly, the dispersion relation for the longitudinal 
waves can be expressed by eliminating Su s , Sn s , and SE Z 
from Eqs. fl64|), ©, and ([SI, of the form 
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(72) 



Like transverse electromagnetic waves it is also true that 
the general relativistic effects enter this dispersion rela- 
tion of longitudinal waves by the ratios of cj/ag. If we 
take the local time dr — a^dt of the FIDO instead of 
a global time t, then the dispersion relation correspond 
to the special relativistic version as it should be accord- 
ing to the Einstein principle of relativity. It is clear that 
the dispersion relation for longitudinal waves is indepen- 
dent of local cyclotron frequency uj cs . If one considers 
the equivalent case to that of SK [24] for an electron- 
positron plasma by taking that the two fluids have the 
same velocity Uq, the same equilibrium density n , and 
are at the same temperature T , one obtains from Eq. 

CI 



1 = 



(uok-^r- v -?(k 



(73) 



In the limit of zero gravity, i.e., when ao — > 1 and uq — > 0, 



the above Eq. ([73]) becomes 

u? = {2oj 2 p + k 2 lg P a / Po ), 



(74) 



which is the SK [24] result. The only difference being 
that 7 g has been set to unity in the SK work [24 . 




Figure 1: Purely real Alfven mode for the electron- 
positron plasma. 



9 Numerical Solutions 

The dispersion relations for the transverse and longi- 
tudinal waves given in Eqs. (|71[) and (|72l) are compli- 
cated enough, even in the simplest cases for the electron- 
positron plasma where both species are assumed to have 
the same equilibrium parameters, and an analytical solu- 
tion is cumbersome and unprofitable. We put the sets of 
two-fluid equations for transverse and longitudinal waves 
separately in the form of a matrix equation in order to 
determine all the physically meaningful modes for each 
waves as follows: 



(A - kI)X = 0. 



(75) 



Where the eigenvalue is chosen to be the wave number k, 
the eigenvector X is given by the relevant set of pertur- 
bations, and I is the identity matrix. The vector A is the 
sum of two matrices An and Aj. The elements of these 
matrices are, respectively, the real and imaginary terms 
in the coefficients of the transverse and longitudinal two- 
fluid equations for perturbation respectively. We need to 
write the perturbation equations in an appropriate form. 
We introduce the following set of dimensionless variables: 
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(76) 



The Su s , Sv s , and u 0s are already dimensionless, but it is 
convenient to defined Su s and Sv s as above consistency. 




Figure 2: Left: one purely imaginary Alfven growth mode 
for the electron-positron plasma. . Right: Another purely 
imaginary Alfven damped mode. 

This is because the choice of input parameters is the same 
for each fluid. By considering transverse electromagnetic 
waves in the gravitational field of Schwarzschild black 
hole, Daniel and Tajima [23] have shown for electron- 
positron plasma that the upper branch represents the 
high frequency electromagnetic waves, which can prop- 
agate in vacuum, but cannot exist for frequencies lower 

then ^J^ujp + u>c, which is the cutoff frequency for this 

plasma and the lower branch represents Alfven waves, 
which cannot exist for a frequency greater then uj c , the 
cyclotron frequency for Alfven waves. Therefore, for 
electron-positron plasma cj* has been chosen as 



Alfven modes, 

high frequency modes 



with ujpi 
and u! r - 



<Jp2 



and us c i = uj C 2, where uj p = ^/0J p iUJp2 
For the case of an electron-ion 



V /W C 1W C 2 

plasma the choice of w* is more complicated matter for 
transverse waves because both the plasma frequency and 
the cyclotron frequency are different for each fluid. So 
it is not clear from the dispersion relation of transverse 
waves what the natural choice of should be. It has 



Figure 3: One purely new real modes for electron-ion 
plasma. Another is shown in Fig. 1. 



been assumed for simplicity that 



where i 



1 (, , 2 

+ 1 



,J c2) 



Alfven modes, 

high frequency modes 



= (2ajp S +w^ s ). The reason for choosing these 
values for is that they reduce to the special relativis- 
tic cutoffs in the zero gravity limits for electron-positron 
plasma. The cutoffs in the special relativistic case are de- 
termined by the dispersion relation in that the solutions 
to the dispersion relation are physical (i.e., Re(k) > 0) 
only for certain frequency regimes. As the dispersion re- 
lation can not be handled analytically, it is difficult to 
determine what the cutoffs should be in the case includ- 
ing gravity. Other similar combination for ui* should not 
make any real difference to the form of the results be- 
cause is really only a scale factor. 

The dimensionless eigenvector for the transverse set 
of equations is 



X: 



t ransverse 



5vi 
5v 2 
SB 
SE 



(77) 



Using Eq. (|76|) . the set of transverse two-fluid equations, 
Eqs. ([65]) . (l66|) . and (|68|) . become the following dimen- 
sionless form: 
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and 
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These equations are now in the required form to be used 
as input to Eq. (|75[) for transverse electromagnetic waves. 
For the longitudinal waves w* and the dimensionless eigen- 
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The dimensionless set of equations for longitudinal waves, 





Figure 5: The left and right imaginary parts show damped 
and growing occur closed to the horizon. The real part of 
these modes is shown in Fig. 1. 

Eqs. ([63]) . (|69| . and dZOj) , can be obtained using Eq. (|76|) 

as 



Figure 4: Top: Real part of the complex conjugate pair 
of Alfven damping and growth modes for the electron-ion 
plasma. Bottom: Left: Imaginary part of the growth 
mode. Right: Imaginary part of the damped mode. 

vector have been taken as 
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and 
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From the above Eqs. (|82|) and ([83]) it is clear that there 
occurs a singularity at the point for which the infall veloc- 
ity equals to the half of the fluid thermal velocity, called 

2 

the transonic radius occur at at so that at = Uq s = , 
for each fluid of species s. The position of the transonic 




Figure 6: The left mode is the real part of high frequency 
mode for the electron-positron plasma. The right mode 
show the plasma is not uniform near the horizon and 
become damped and growing. 

radius for each fluid is principally depend on the fluid 
temperature, the limiting temperature of each fluid at 
the horizon, which determines the temperature at any 
given radius. Therefore, the transonic radius plays a sig- 
nificant role for the longitudinal waves and its influence 
on the on the waves will be clearly shown by some of 
the longitudinal waves modes. Equations (|82|) . (|83|) . and 
(|84|) are the required equations to be used as input to 
Eq. (|75|) for longitudinal waves which can be put as 
(A — kI)X = 0. The eigenvalues k of the complex ma- 
trix A have been calculated to draw the modes for the 



transverse electromagnetic and longitudinal waves using 
MATLAB. 
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Figure 7: The left mode is the real part of the two complex 
conjugate high frequency mode for the electron-positron 
plasma. The right mode show the plasma is not uniform 
near the horizon and become growth near the horizon. 



9.1 Results 

We have considered both the electron-positron plasma 
and the electron-ion plasma. The limiting horizon values 
for the electron-positron plasma are taken to be 

n Hs = 10 18 cm- 3 , T Hs = 10 10 K,B# = 3x 10 6 G, and lg 

For the electron-ion plasma, the ions are essentially non- 
relativistic, and the limiting horizon values are chosen to 
be 

nm = 10 18 cnT 3 , T H1 = 10 10 K,n ff2 = 10 15 cm" 3 , 

and T H2 = 10 12 K. 

The equilibrium magnetic field has the same value as it 
has for the electron-positron case. The limiting horizon 
temperature for each species has been chosen as derived 
by Colpi et al. [36] from studies of two temperature 
models of spherical accretion onto black holes. Those 
values of the limiting horizon densities and the limiting 
horizon field are arbitrarily chosen which appear to be 




Figure 8: Both the modes are stable for all the values 
of ao but very closed to the horizon they are unstable 
and the damped and growing occur for all the frequencies. 
The real part is shown in Fig. 6. 

not inconsistent with current ideas. The gas constant 
has been chosen as "f g = |. 

9.2 Alfven Modes 

9.2.1 Electron-positron Plasma 

For the special relativistic electron-positron plasma, only 
one purely real Alfven mode was found to exist by SK 
[24] . Because both the left and right circularly polarized 
modes were described by the same dispersion relation. 
Our work presented in this paper also shows one purely 
real Alfven mode with two new purely imaginary con- 
jugate modes to exist for the same plasma due to the 
gravitational field of Schwarzschild black hole. The real 
mode shown in Fig. 1 is a common mode for special and 
general relativistic electron-positron plasma but the two 
purely imaginary modes shown in Fig. 2 are respectively 
damped and growing and the damped and growing rate 
are clearly frequency independent, depend only on lapse 
function ao . The damped mode demonstrates the energy 
drain from the waves by the gravitational field and the 
growth mode indicates that the gravitational field is, in 
fact, feeding energy into the waves. Both the modes are 



stable for all the frequencies and at all radial distances 
from the event horizon through a^. Since we are using 
the convention e ^ k ^ dz = e 1 m<k)+iim(k)]dz j the damp- 
ing corresponds to Im(k) > and growth to Im(k) < 0. 
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Figure 9: Instability occur very closed to the horizon for 
electron-ion plasma. The real part is shown in Fig. 7. 



9.2.2 Electron-ion Plasma 

In the case four modes are found to exist, two of which 
are purely real and others are complex conjugate modes. 
The mode shown in Fig. 1 for electron-positron plasma 
also found to exist for electron-ion plasma. The other 
real mode shown in Fig. 3 is a new mode for electron-ion 
plasma. Thus, for electron-ion plasma two purely real 
modes exist due to the difference in mass and density 
factors as between the positrons and ions. These two 
modes coalescence with a single mode for the special rel- 
ativistic electron-positron plasma as investigated by SK 
[23] . The other two modes shown in Fig. 4, are a complex 
conjugate pair and are significantly damped and growing, 
respectively. The imaginary parts of this complex conju- 
gate modes are equivalent to the two new purely imagi- 
nary modes we have found for electron-positron case dis- 
cussed above with larger damping and growth rates. The 
differences in the magnitudes of the cyclotron frequen- 
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Figure 10: Purely real longitudinal mode for the electron- 
positron plasma. 

ties u) c i and u) C 2 for the last two modes apparently lead 
to take the frequencies from their negative (and there- 
fore unphysical) values for the electron-positron case to 
positive physical values for the electron-ion case. It is 
evident that these two modes are also stable and the 
damping and growth rates are independent of frequency, 
dependent only on the distance from the black hole hori- 
zon. 

9.3 High Frequency Transverse Modes 

9.3.1 Electron-positron Plasma 

Four high frequency electromagnetic modes are found to 
exist for the electron-positron plasma. Two modes shown 
in Fig. 5 have same real parts as shown in Fig. 1. The ef- 
fect of the general relativistic term is clearly shown here. 
For lower frequencies the plasma are unstable as moved 
toward the horizon and the damping and growth occurs 
as a — > 0. It then appears, at a distance from the horizon 
corresponding to a < 0.4 and Co < 6 that energy is no 
longer fed into wave mode by the gravitational field but 
begins to be drained from the waves. The third mode, 
shown in Fig. 6, is also unstable near the event horizon. 
For lower frequencies the damping and growth rates are 
very high, but for higher frequencies the damped and 
growing rates are smaller. The fourth mode, shown in 
Fig. 7 is stable for a > 0.6 and for higher frequencies, 
but unstable for a < 0.6 and for lower frequencies. That 
is to say that, the damped and growing of these modes 
are depended on both the wave frequency and lapse func- 
tion a . 

9.3.2 Electron-ion Plasma 

Like the electron-positron plasma, the electron-ion plasma 
admits four high frequency modes. These are illustrated 
in Figs. 8 and 9. The Fig. 8 shows two complex modes 
having real parts equal and are similar with the modes 
shown in Fig. 1 for electron-positron plasma. Both 
the imaginary part of these two modes are stable for all 



frequencies and at all distances from the horizon corre- 
sponding to ao > 0.1, but the solution becomes unstable 
for ao < 0.1. This is also true for the remaining two 
modes shown in Fig. 9. The real part of this two modes 
are similar with the two real part of the modes shown in 
Figs. 6 and 7, respectively. Very closed to the horizon the 
growth and decay rates depend on both frequency and 
lapse function. The damping and growth are not clearly 
evident like the corresponding high frequency modes for 
electron-positron plasma. 




Figure 11: Left: Real part of a longitudinal complex 
mode for the electron-positron plasma. Right: The 
imaginary part shows the interplay between frequency Co 
and a near the singularity and near the event horizon. 

9.4 Longitudinal Modes 

9.4.1 Electron-positron Plasma 

The longitudinal modes are not split into high and low 
frequency domain as was done for transverse electromag- 
netic waves because the same modes are exist in the low 
frequency, < uj < 1, and high frequency, 1 < Co < 10, 
domains. So, Co has been chosen here as < Co < 10. For 
electron-positron plasma, there exist five modes one of 
which is purely real mode shown in Fig. 10. This is also 
true for the special relativistic case investigated by SK 
where only one high frequency mode was found to exist 
for electron-positron plasma. For this mode transonic ra- 
dius occurs at about a t ~ 0.82. The second mode shown 




Figure 12: Left: Real part of a longitudinal complex 
mode for the electron-positron plasma. Right: The 
imaginary part clearly shows the interplay between fre- 
quency Co and ao near the transonic radius, at ~ 0.82, 
and near the event horizon, ao — > 0. 

in Fig. 11 is a complex mode with real part equal to the 
real modes shown above figure. In the imaginary part of 
this mode some interplay between frequency uj and ao 
occur near the singularity and near the event horizon. 
The third modes shown in Fig. 12 also shows some in- 
terplay between u) and lapse function. The fourth mode 
shown in Fig. 13 is growth for ao < at and damped 
for ao > at- This mean that energy is drained from the 
wave rather than being fed into it by the gravitational 
field. The fifth mode is almost the opposite of the pre- 
vious mode in that it is growth mode for ao < at and 
damped for ao > at. 

9.4.2 Electron-ion Plasma 

As for electron-positron case, here five low frequency 
modes exist three of which are purely real and others 
are complex. Since the two-fluid (electron and ion) tem- 
peratures near the horizon are different, their transonic 
radii are different and are occurred at a t i ~ 0.82 and 
a t 2 ~ 0.99 respectively. Figure 14 shows one purely real 
mode for electron-positron plasma. The other two purely 
real modes are shown previously in Fig. 10 and the top 
of the Fig. 12. The fourth mode shown in Fig. 15, shows 



both damping and growth. It is stable for ao < an and 
growth for an < ao < a*2 and then, for ao > an it 
transition back to damped. The fifth mode, not shown 
here, is almost the opposite of the fourth mode and is 
stable for ao < an and damped for an < &o < &t2 and 
then, for ao > ati it transition back to growth. These 
two modes also clearly show the influence of the tran- 
sonic radii and again the fact that energy is being fed 
into the wave between the transonic radii but is drained 
from the wave very closed to the horizon. 
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Figure 13: Left: Real part of a longitudinal complex 
mode for the electron-positron plasma. Right: Imagi- 
nary part shows damped corresponding to ao < at. 



10 Concluding Remarks 

The prime concern of this study has been exclusively the 
investigation, within the WKB approximation, of trans- 
verse electromagnetic and longitudinal waves in a two- 
fluid plasma closed to the event horizon of Schwarzschild 
black hole. We have derive the local dispersion relations 
for such types of waves and shown that the general rel- 
ativistic effects enter these dispersion relations only in 
the ratio w/ao and with dr = a$dt the transverse and 
longitudinal waves dispersion relations reduced to the 
special relativistic version as it should be according to 
the Einstein relativity principle. We solve the dispersion 
relations numerically for the wave number k. For the 



accordance with the results of analytical studies carried 
out respectively by Mikhailovskii [33] and Rahman [31] . 








Figure 14: Three purely real longitudinal modes exist for 
electron-ion plasma one of which is shown above and the 
another two are shown previously in Fig. 10 and the left 
of the Fig. 12. 

electron-positron plasma, the damping and growth rates 
are smaller in general, by several orders of magnitude, 
compared with the real components of the wave number. 
On the contrary, for the electron-ion plasma, modes with 
significant damping and growth rates are found to exist, 
in particular, for the Alfven waves. Obviously, the damp- 
ing and growth rates are frequency independent for the 
Alfven waves, but they are solely dependent on the ra- 
dial distance from the horizon through the lapse function 
olq. However, this is not the case for the high frequency 
waves. The rate of damping or growth in that case is de- 
pendent on both frequency and radial distance from the 
horizon. These results are essentially agree with the re- 
sults of SK [24] for ultrarelativistic limit and of Rahman 
[31 and Daniel and Tajima 23 for general relativity. For 
very low or negligible frequency, the plasmas are unsta- 
ble, i.e., damped and growing modes exist as mansion 
above. Similar instabilities of a relativistic plasma were 
found by Mikhailovskii 33 , and Zaslavskii and Moiseev 
[?] . The same conclusion also follows from the numerical 
calculations of Buzzi et al. [25] [26] . 

For longitudinal waves the influence of transonic ra- 
dius for each of the fluid species has been clearly ev- 
ident. Like transverse electromagnetic waves damping 
and growth modes are found for each fluid species. The 
damped and growing rates for some of the modes depen- 
dent on frequency but for all the modes, the damping and 
growth rates are dependent on the radial distance from 
the horizon through the lapse function cto . For each fluid 
species, all the real transverse electromagnetic modes in 
the ao domain are equivalent to all the real longitudinal 
modes in the ao < &t domain and give the same results 
for the special relativistic case investigated by SK [2~i] . 
The presence of damped modes demonstrates that, at 
least in this approximation, energy is being drained from 
the waves by the gravitational field. Since the majority 
of the modes show growth rates, the gravitational field is, 
in fact, feeding energy into the waves. This results are in 




Figure 15: Left: Real part of a longitudinal complex 
mode for the electron-ion plasma. Right: Imaginary 
part of longitudinal mode which is stable for ao < at 
and which shows growing and damped corresponding to 
an < ao < a t 2 and ao > a t 2, respectively. 
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